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7.4 RESOLUTION PERFORMANCE ASSESSMENT FOR 
QUADRATIC TFDs 0 

7.4.1 Selecting and Comparing TFDs 

Quadratic time-frequency distributions (TFDs) are effective tools for extracting in- 
formation from a non-stationary signal, such as the number of components, their 
durations and bandwidths, components’ relative amplitudes and instantaneous fre- 
quency (IF) laws (see Chapters 1 and 2). The performance of TFDs depends on 
the type of signal (see Chapter 3) [1,2]. For example, in the case of a monocompo- 
nent linear FM signal, the Wigner-Ville distribution is known to be optimal in the 
sense that it achieves the best energy concentration around the signal IF law (see 
Article 2.1 for more details) [1]. 

In applications involving multicomponent signals, choosing the right TFD to 
analyze the signals is an immediate critical task for the signal analyst. How best to 
make this assessment, using current knowledge, is the subject of this article. 

Let us, for example, consider a multicomponent whale signal, represented in the 
time-frequency domain using the Wigner-Ville distribution, the spectrogram, the 
Choi-Williams distribution, the Born-Jordan distribution, the Zhao- Atlas-Marks 
(ZAM) distribution, and the recently introduced B-distribution [3] (see Fig. 7.4.1). 

To determine which of the TFDs in Fig. 7.4.1 “best” represents this whale signal 
(i.e. which one gives the best components’ energy concentration and best interfer- 
ence terms suppression, and allows the best estimation of the components’ IF laws) 
one could visually compare the six plots and choose the most appealing. The spec- 
trogram and the B-distribution, being almost free from the cross-terms, seem to 
perform best. 

The performance comparison based on the visual inspection of the plots becomes 
more difficult and unreliable, however, when the signal components are closely- 
spaced in the time-frequency plane. To objectively compare the plots in Fig. 7.4.1 
requires to use a quantitative performance measure for TFDs. There have been 
several attempts to define objective measures of “complexity” for TFDs (see Sec- 
tion 7.3.1). One of these measures, the Renyi entropy given in [4], has been used by 
several authors in preference to e.g. the bandwidth-duration product given in [1], 
The performance measure described in this article, unlike the Renyi entropy, is a 
local measure of the TFD resolution performance, and is thus more suited to the 
selection problem illustrated by Fig. 7.4.1. This measure takes into account the 
characteristics of TFDs that influence their resolution, such as energy concentra- 
tion, components separation, and interference terms minimization. Methodologies 
for choosing a TFD which best suits a given signal can then be developed by optimiz- 
ing the resolution performance of considered TFDs and modifying their parameters 
to better match application-specific requirements. 
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(b.boashash@qut.edu.au, v.sucic@qut.edu.au). Reviewers: W. J. Williams and LJ. Stankovic. 
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(a) Wigner-Ville 




(b) Spectrogram (Han- 
ning, L = 511) 





(d) Born-Jordan (e) ZAM (a = 2) (f) B (/3 = 0.01) 

Fig. 7.4.1: TFDs of a multicomponent whale signal. 



7.4.2 Performance Criteria for TFDs 

In the case of monocomponent FM signals, the best TFD is that which maximizes 
energy concentration about the signal instantaneous frequency. This is achieved by 
minimizing component sidelobe amplitude A s relative to mainlobe amplitude A m , 
and mainlobe bandwidth B relative to central frequency / (see Fig. 7.4.2). 

The instantaneous concentration performance of a TFD may thus be quantified 
by the measure p expressed as: 



P(t) 



A s {t) Bjt) 

A m{t) fit) 



(7.4.1) 



A good performance is characterized by a small value of the measure p. For ex- 
ample, for the Wigner-Ville distribution of a linear FM signal with infinite duration, 
the bandwidth B and the sidelobe amplitude A s are zero [1], and we obtain p = 0. 

For multicomponent FM signals, the performance of a TFD can be quantitatively 
assessed in terms of: 



• the energy concentration of the TFD about the respective instantaneous fre- 
quency of each component, as expressed by Eq. (7.4.1), and 

• the components resolution, as measured by the frequency separation of the 
components’ mainlobes, including the effect of cross-terms. 
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Fig. 7.4.2: Slice of a TFD of a monocomponent signal at time t = to. The dominant peak is the 
component, while the other peaks are the sidelobes. For clarity of presentation, we limit ourselves to 
measuring the mainlobe bandwidth at 0.71 of the component normalized amplitude Am. 




Fig. 7.4.3: Diagram illustrating the resolution of a two-component signal in the absence of cross-terms. 
The lobes are clearly distinguished from each other; the components are said to be resolved. 



For stationary signals, the frequency resolution in a power spectral density es- 
timate of a signal composed of two single tones, /i and / 2 (see Fig. 7.4.3), may be 
defined as the minimum difference f 2 — /i for which the following inequality holds: 

fi + £i/2 <h- B 2 /2 , h < f 2 (7-4.2) 

where B\ and B 2 are the respective bandwidths of the first and the second sinusoid. 

In the case of non-stationary signals, for a TFD p z (t, f) of a two-component 
signal, the above definition of resolution is valid for every time slice of a cross-term- 
free TFD, such as the spectrogram. However, for TFDs exhibiting cross-terms, we 
need to take into account the effect of cross-terms on resolution. 

A slice of a typical quadratic TFD, with components clearly resolved, is shown 
in Fig. 7.4.4(a), where Bi(to), fi(to), A Sl (to) and A mi (to) represent respectively 
the instantaneous bandwidth, the IF, the sidelobe amplitude and the mainlobe 
amplitude of the first component at time t = t$. Similarly, B 2 (to), / 2 (to), A S2 (to) 
and A m2 (t 0 ) represent the instantaneous bandwidth, the IF, the sidelobe amplitude 
and the mainlobe amplitude of the second component at the same time to- The 
amplitude A x (to) is that of the cross-term. An example of a quadratic TFD with 
non-resolved components is shown in Fig. 7.4.4(b). 
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(a) The two dominant peaks are the (b) The two components and the 
(resolved) signal components, the cross-term have merged into a single 

middle peak is the cross-term, and peak; we say that the components 

the other peaks are the sidelobes are not resolved 

Fig. 7.4.4: Slice of a TFD of a two-component signal taken at time t = to- 



7.4.3 Resolution Performance Measure for TFDs 

Eq. (7.4.2) and Fig. 7.4.4(a) suggest that the resolution performance of a TFD for a 
neighboring pair of components in a multicomponent signal may be defined by the 
minimum difference D(t) = / 2 (f) — /i(t) for which we still have a positive separation 
measure S(t) between the components’ mainlobes centered about their respective 
IFs, fi(t) and / 2 (f). For best resolution performance of TFDs, 5(f) should be as 
close as possible to the true difference between the actual IFs. 

The components’ separation measure 5(f) is expressed as [5]: 








(7.4.3) 



The resolution also depends on the following set of variables, all of which should 
be as small as possible: 

(a) the normalized instantaneous bandwidth of the signal component Bk(t)/ fk{t), 
k = 1, 2, which is accounted for in 5(f) (Eq. (7.4.3)), 

(b) the ratio of the sidelobe amplitude |H Sfe (f)| to the mainlobe amplitude \A mk (t)\, 
/c = 1,2, of the components, and 



(c) the ratio of the cross-term amplitude |A„(t)| to the mainlobe amplitudes of 
the signal components \A mk (t)\, k = 1,2. 

It follows that the best TFD for multicomponent signals analysis is the one that 
concurrently minimizes the positive quantities (a), (b), (c), and maximizes 5(f). 

Hence, by combining the above variables, expressions for a measure P(f) of the 
resolution performance of a given TFD can be defined. Two have been proposed 
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in [5], among these a normalized performance measure expressed as: 



!/ 


A s (t) 


1 


A x {t) 


3 l 


A m (t) 


+ 2 


A m {t) 




(7.4.4) 



where, for a pair of signal components, A m ( t ) and A s (t) are respectively the average 
amplitudes of the components’ mainlobes and sidelobes, A x (t) is the cross-term 
amplitude, S(t), defined by Eq. (7.4.3), is a measure of the components’ separation in 
frequency, and D(t ) = / 2 (f) — f\ (t) is the difference between the components’ actual 
IFs. The algorithm presented in [6] describes how the parameters in Eq. (7.4.4) are 
measured in practice. 

The measure P(t) is close to 1 for well-performing TFDs and 0 for poorly- 
performing ones. Therefore, when designing a TFD we want to maximize P(t) in 
order to reduce the cross-terms, while preserving the components’ resolution. 

In some applications involving real-life signals, we may need to better discrimi- 
nate between different TFDs resolution performances in a given set of K TFDs. In 
this case, a suitable alternative to P(t) that was proposed in [7] could be used. It 
is expressed as: 



Mj(t) = 1-- 





As j ( t ) 
Am j (t) 




max | 

l<k<K 


; 


A s k (t) 

Am k (t) 


) 



A X j (t) 

A^Aij 



Dj(t) 



max 

l<k<K 



( A Xk ( t ) \ 

\ Am k (t) ) 



( B k (t) 

max 

l<k<K \ Dk ^ t ) 



(7.4.5) 

where M 7 (i) (1 < j < K) is the resolution performance measure of the j-th TFD, 
and B is the average instantaneous bandwidth of the components mainlobes. The 
measure M(t ) is used in Section 7.4.5 to compare the performances of quadratic 
TFDs of a real-life signal, as it discriminates better than the measure P(t) for 
real-life signals [8] . 



7.4.4 Application to the Selection of the Optimal TFD 
for a Given Multicomponent Signal 

A methodology for selecting the optimal TFD for resolving closely-spaced compo- 
nents in a multicomponent signal involves then the following steps: 

1. Define a set of comparison criteria describing the information sought from 
TFDs (Section 7.4.2). 

2. Objectively measure the resolution performance of TFDs based on these cri- 
teria (use the measure P defined by Eq. (7.4.4)). 

3. Optimize each TFD to match the criteria as close as possible [5,6]: Select as 
the optimal TFD kernel parameter value the one which maximizes the over- 
all performance measure Poverain taken to be the mean of the instantaneous 
measures P in a time interval of interest. 
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Table 7.4.1: Optimization results for the TFDs of signal s(t) defined by Eq. (7.4.6). 



TFD 


Optimal value of the kernel parameter 


P Dverall 


Spectrogram 


Bartlett window, length 31 


0.86 


Wigner-Ville 


N /A in this case 


0.62 


Choi-Williams 


(7=1 


0.82 


Born-Jordan 


N /A in this case 


0.81 


Zhao-Atlas-Marks (ZAM) 


a = 2 


0.67 


Modified B 


/? = 0.04 


0.88 



4. Quantitatively compare TFDs and select the best one: An optimized TFD 
which has the largest value of P 0 veraii is selected as best for representing the 
given signal in the joint time-frequency domain. 

Example: We define the following two-component signal in noise: 

s(t) = si(t) + s 2 (t) + n{t) 

= cos (27 r (0.1 1 T f t 2 )) + cos (27r (0.2 1 + f t 2 )) + n{t) (7.4.6) 

where a = 0.0016 is the component bandwidth-duration ratio (duration T = 128), 
and n(t) is additive white Gaussian noise with signal-to- noise ratio SNR = 10 dB. 
The sampling frequency is f s = 1 Hz. 

The signal s(t) is analyzed in the time-frequency domain using the following 
TFDs: the spectrogram, the Wigner-Ville distribution, the Choi-Williams distribu- 
tion, the Born-Jordan distribution, the Zhao-Atlas-Marks (ZAM) distribution, and 
the Modified B-distribution [9]. 

To find the optimal TFD for resolving the two components of s(f), we first 
find the optimal values of the TFDs kernel parameters, as described in the above 
methodology. The Wigner-Ville distribution and the Born-Jordan distribution have 
no “smoothing” parameters, hence do not need optimizing. The optimized TFD 
with the largest P 0 veraii among the considered TFDs is then selected as optimal 
for representing s(t). Table 7.4.1 lists the results of the optimization process, and 
it shows that the signal optimal TFD is the Modified B-distribution with the pa- 
rameter (3 = 0.04. All optimized TFDs are plotted in Fig. 7.4.5. From the signal 
optimal TFD important signal parameters can be measured (see Table 7.4.2). In 
addition, by optimizing components’ concentration and resolution, more accurate 
components IF laws’ estimates are obtained from the peaks of the optimal TFD’s 
dominant ridges in the time-frequency plane [1] (see Fig. 7.4.6). 

7.4.5 Use of the Performance Measure in Real-Life Situations 

The methodology defined in this section enables to select a real-life signal best- 
performing TFD in an objective, automatic way. Its use should make time- frequency 
techniques more applicable in practice (e.g. machine condition monitoring described 
in Articles 15.2 and 15.6, or other applications presented in Part V). 
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(a) Spectrogram (Bartlett, 
L = 31) 




(b) Wigner-Ville 






(e) ZAM (a = 2) 




(f) Modified B (0 = 0.04) 



Fig. 7.4.5: Optimized TFDs of signal s(t ) defined by Eq. (7.4.6). 



Table 7.4.2: Parameters of si(t) and S 2 (f) (Eq. (7.4.6)) measured from the signal s(t) optimal TFD 
(Modified B-distribution, 6 = 0.04). The values shown are the averages over t £ [32,96]. 



Parameter 


Component si(t) 


Component S 2 (t) 


Instantaneous bandwidth B 


0.0194 


0.0195 


Mainlobe amplitude \A m \ 


1.0002 


0.9574 


Sidelobe amplitude A s 


0.0900 


0.0858 


Cross-term amplitude \A X \ 


0.1503 



The methodology consists of the following steps: 

1. Represent the signal in the time- frequency domain with a quadratic TFD, i.e. a 
smoothed Wigner-Ville distribution (see Article 3.2). Following the approach 
described in Article 5.7, we smooth the WVD in both time t and lag r with 
the Hanning window of length equal to a quarter of the signal duration. This 
time-frequency smoothing is intended to suppress the WVD inner and outer 
artifacts, while preserving components time-frequency features (see Article 4.2 
for more details). 

2. For the different time instants of the smoothed WVD, select the two closest 
dominant peaks in the frequency direction. To achieve the best resolution of 
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(a) Modified B (/3 = 0.04) (b) Spectrogram (Bartlett, L = 31) 

Fig. 7.4.6: Comparison of the measured (dashed) and true (solid) IF laws of the component Sl(t) 
(left) and S 2 (t ) (right) of the signal s(t ) defined by Eq. (7.4.6). The mean-square-errors (MSEs) of the 
IF estimates obtained from the peaks of the signal optimized Modified B-distribution (best-performing 
TFD) are given in (a), and those obtained from the peaks of the signal optimized spectrogram (second 
best TFD) in (b). 



the signal components, the best resolution of the two closest components at an 
observed time instant is sufficient [10]. Note that if a signal is monocomponent 
or no components exist at a particular time, this time instant is not considered. 

3. For the selected pairs of components, optimize different TFDs using the reso- 
lution performance measure M defined by Eq. (7.4.5). The measure M is used 
over P since it is a better discriminator of real-life signals TFDs resolution 
performances [11]. The kernel parameter value, which from a set of different 
values considered, maximizes the overall performance measure (the mean of 
M over the observed times) is selected as the kernel parameter optimal value. 

4. Calculate the measure M of the optimized TFDs for each of the selected pairs 
of signal components. The TFD which maximizes the average (over time) M 
is selected as the signal best-performing TFD among the considered TFDs. 

Example: To illustrate how to use this methodology in practice, let us find the best- 
performing TFD for the Noisy Miner ( Manorina melanocephala) song signal. The 
same TFDs we considered in the synthetic signal optimal TFD selection example 
will be considered in this real-life signal example. 

We start by representing the signal in the time-frequency domain with the 
Wigner-Ville distribution smoothed in time and lag with the Hanning windows 
of length L = 3501 (Fig. 7.4.7(a)). 

For each time instant of the smoothed WVD we then identify the pair of closest 
components. From Fig. 7.4.7(b) we can see that different components form such 
pairs at different times. 

Next, as described in steps 3 and 4 of the above-defined methodology, the six 
considered TFDs are first optimized, after which their resolution performances are 
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(a) Smoothed WVD (Han- (b) Pairs of closest compo- 
ning, L = 3501) nents 




(c) Modified B (/? = 0.91) 



Fig. 7.4.7: Optimization of the bird song signal TFDs. The signal Smoothed WVD is shown in (a), 
and the pairs of its closest components in (b), with the overall closest pair (at time t = 4295) marked 
by the squares. The signal optimized Modified B-distribution is shown in (c). 



Table 7.4.3: Optimization and comparison results for the TFDs of the Noisy Miner song signal. The 
values of M overa n indicate that the spectrogram performs better than most traditional TFDs in this 
case. Only the Modified B-distribution performs better than the spectrogram and all others. 



TFD 


Optimal value of the kernel parameter 


-^overall 


Spectrogram 


Bartlett window, length 3501 


0.90 


Wigner-Ville 


N/A in this case 


0.50 


Choi-Williams 


cr = 0.004 


0.74 


Born-Jordan 


N/A in this case 


0.65 


Zhao-Atlas-Marks (ZAM) 


a = 2 


0.63 


Modified B 


/3 = 0.91 


0.93 



evaluated using the measure M . Table 7.4.3 shows the signal TFDs kernel pa- 
rameters optimization and the TFDs resolution performance results. The Modified 
B-distribution for /3 = 0.91, plotted in Fig. 7.4.7(c), is found to have the largest 
value of M OV e r aii (the mean of M over the time instants). Therefore, we select this 
TFD as best to represent the Noisy Miner song signal in the time-frequency plane. 



7.4.6 Summary and Conclusions 

This article defines a measure for assessing the resolution performance of quadratic 
TFDs in separating closely-spaced components in the time-frequency domain. The 
measure takes into account key attributes of TFDs, such as components’ mainlobes 
and sidelobes, and cross-terms. The introduction of this measure allows to quantify 
the quality of TFDs instead of relying solely on visual inspection of plots. The 
resolution performance measure also allows for selecting the optimal TFD in a given 
practical application, and improving methodologies for designing high resolution 
quadratic TFDs, such as the Modified B-distribution. 
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